arXiv:1503.02996vl [hep-ph] 10 Mar 2015 


Higgs Decays into Charginos and Neutralinos 
in the Complex MSSM: A Full One-Loop Analysis 

S. HeinemeyerIiI and C. Schappacher® 

^Instituto de Fisica de Cantabria (CSIC-UC), Santander, Spain 

^Institut fiir Theoretische Physik, Karlsruhe Institute of Technology, 
D-76128 Karlsruhe, Germany 


Abstract 

For the search for additional Higgs bosons in the Minimal Supersymmetric Standard 
Model (MSSM) as well as for future precision analyses in the Higgs sector a precise 
knowledge of their decay properties is mandatory. We evaluate all two-body decay 
modes of the Higgs bosons into charginos and neutralinos in the MSSM with complex 
parameters (cMSSM). The evaluation is based on a full one-loop calculation of all decay 
channels, also including hard QED radiation. We restricted ourselves to a version 
of our renormalization scheme which is valid for |Mi| < |M 2 |,|r| and M 2 7 ^ r to 
simplify the analysis, even though we are able to switch to other parameter regions. 
The dependence of the Higgs boson predictions on the relevant cMSSM parameters 
is analyzed numerically. We find sizable contributions to many partial decay widths. 
They are roughly of 10% of the tree-level results, but can go up to 20% or higher. The 
full one-loop contributions are important for the correct interpretation of heavy Higgs 
boson search results at the LHC and, if kinematically allowed, at a future linear e'’'e“ 
collider. It is planned to implement the evaluation of the branching ratios of the Higgs 
bosons into the Fortran code FeynHiggs, together with an automated choice of the 
renormalization scheme valid for the full cMSSM parameter space. 
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1 Introduction 


One of the most important tasks at the LHC is to search for physics effects beyond the 
Standard Model (SM), where the Minimal Supersymmetric Standard Model (MSSM) [TH3] 
is one of the leading candidates. Supersymmetry (SUSY) predicts two scalar partners for 
all SM fermions as well as fermionic partners to all SM bosons. Another important task 
of the LHC is the investigation of the mechanism of electroweak symmetry breaking. The 
most frequently studied realizations are the Higgs mechanism within the SM and within the 
MSSM. Contrary to the case of the SM, in the MSSM two Higgs doublets are required. This 
results in hve physical Higgs bosons instead of the single Higgs boson in the SM. In lowest 
order these are the light and heavy CP-even Higgs bosons, h and H, the CP-odd Higgs 
boson, A, and two charged Higgs bosons, Within the MSSM with complex parameters 
(cMSSM), taking higher-order corrections into account, the three neutral Higgs bosons mix 
and result in the states p {i = 1, 2, 3) [1H7]. The Higgs sector of the cMSSM is described at 
the tree-level by two parameters: the mass of the charged Higgs boson, Mh± , and the ratio 
of the two vacuum expectation values, tan/3 = tp = V 2 /V 1 . Often the lightest Higgs boson, 
hi is identihed with the particle discovered at the LHC [SllH] with a mass around ~ 125 GeV. 
If the mass of the charged Higgs boson is assumed to be larger than ~ 200 GeV the four 
additional Higgs bosons are roughly mass degenerate, Mh± k, nth^ k, and referred to as 
the “heavy Higgs bosons”. Discovering one or more of those additional Higgs bosons would 
be an unambiguous sign of physics beyond the SM and could yield important information 
about their possible supersymmetric origin. 

If SUSY is realized in nature and the charged Higgs-boson mass is Mh± 1-5 TeV, 
then the heavy Higgs bosons could be detectable at the LHC [lOllTT] (including its high 
luminosity upgrade, HL-LHC, see Ref. [12] and references therein) and/or at a future linear 
e’''e“ collider such as the ILC [T31 - IT5] or CLIC [16]. (Results on the combination of LHC and 
ILC results can be found in Ref. El-) The discovery potential at the HL-LHC goes up to 
0{\ TeV) for large tan/3 values and somewhat lower at low tan/3 values. At an e^e“ linear 
collider the heavy Higgs bosons are pair produced, and the reach is limited by the center 
of mass energy, Mh± < \fs/2, roughly independent of tan/3. Details about the discovery 
process(es) depend strongly on the cMSSM parameters (and will not be further discussed in 
this paper). 

In the case of a discovery of additional Higgs bosons a subsequent precision measure¬ 
ment of their properties will be crucial determine their nature and the underlying (SUSY) 
parameters. In order to yield a sufficient accuracy, one-loop corrections to the various Higgs- 
boson decay modes have to be considered. Decays to SM fermions have been evaluated at 
the full one-loop level in the cMSSM in Ref. [18], see also Ref. [19] as well as Refs. [2011^ 
for higher-order SUSY corrections. Decays to (lighter) Higgs bosons have been evaluated 
at the full one-loop level in the cMSSM in Ref. [H], see also Refs. [22l|23]. Decays to SM 
gauge bosons (see also Ref. [21]) can be evaluated to a very high precision using the full SM 
one-loop result [25] combined with the appropriate effective couplings [26] . The full one-loop 
corrections in the cMSSM listed here together with resummed SUSY corrections have been 
implemented into the code FeynHiggs [261130] . Corrections at and beyond the one-loop level 
in the MSSM with real parameters (rMSSM) are implemented into the code HDECAY [3T1I32] . 
Both codes were combined by the LHC Higgs Cross Section Working Group to obtain the 
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most precise evaluation for rMSSM Higgs boson decays to SM particles and decays to lighter 
Higgs bosons [5B] . 

The (heavy) MSSM Higgs bosons can (if kinematically allowed) also decay to SUSY 
particles, i.e. to charginos, neutralinos and scalar fermions. In Ref. [M] it was demonstrated 
that the SUSY particle modes can dominate the decay of the heavy Higgs bosons. The 
lightest neutral Higgs boson, on the other hand, can have a substantial branching ratio into 
the lightest neutralino, hi —)■ XiXi, where the Xi constitutes the Dark Matter candidate in 
the MSSM |32]. Bounds on m^o often assume an underlying SUSY grand unihed theory, 
based on a simple Lie group. Dropping these assumptions hardly any bound on m^o can 
be placed directly (see, e.g.. Ref. [3B] and references therein), and the decay hi —)> XiXi 
is kinematically possible. In order to determine the Dark Matter properties a precision 
measurement of this process at the LHC or a future e’''e“ collider will be necessary. 

Higher-order contributions to MSSM Higgs boson decays to scalar fermions have been 
evaluated in various analyses over the last decade. For calculations in the rMSSM, see 
Refs. [371139] and references therein. More recently, the results of Ref. [37| were made public 
in the code HFOLD [10] , using a pure DR renormalization for the calculation. In Ref. PH the 
0(as) corrections to Higgs boson decays to scalar quarks were re-analyzed and included into 
the code HDECAY. Within the cMSSM a full one-loop calculation of Higgs boson decays to 
scalar fermions has recently been published in Ref. [12| and will be included into the code 
FeynHiggs. These results were obtained in a renormalization scheme [2611131118] . which has 
been shown to yield stable results over nearly the full cMSSM parameters space. In this 
work we take another step in the direction of completion of the calculation of all two-body 
decays at the one-loop level in the cMSSM in this stable and reliable renormalization scheme: 
we calculate ah two-body decay modes of the Higgs bosons to charginos and neutralinos in 
the cMSSM. More specihcahy, we calculate the full one-loop corrections to the partial decay 
widths 


r(/ii ^ Xc xt') 

{i = 1,2,3; c,c' = 1,2), 

( 1 ) 

-X xlxl) 

{i = 1, 2, 3; n, n' = 1, 2, 3,4), 

( 2 ) 


(n = 1,2, 3,4; c = 1,2), 

(3) 


where xf iXn) denotes the charginos (neutralinos). While we have calculated the decay of 
ah Higgs bosons, in the numerical evaluation below, we will concentrate on the heavy Higgs 
bosons, h 2,3 and but also show results for hi —)■ X°iXv 

The evaluation of the channels Eqs. o - o is based on a full one-loop calculation, 
i.e. including electroweak (EW) corrections, as well as soft and hard QED radiation. For 
“mixed” decay modes, we evaluate in addition the two “CP-versions” (c ^ c') of Eq. ([T]) 
and the two “CP-versions” of Eq. (|3|), which give different results for non-zero complex 
phases. We restricted ourselves to a version of our renormalization scheme which is valid for 
\Mi\ < IM 2 I, |/i| and M 2 7 ^ n (where Mi and M 2 denote the soft SUSY-breaking parameter 
of the U{1) and SU{2) gauginos, and p is the Higgs mixing parameter) to simplify the 
analysis, even though we are able to switch to other parameter regions, see the discussion in 
Refs. [^GtiiS] (see also Ref. |19]). 

Higher-order contributions to MSSM Higgs boson decays to charginos and neutralinos 
have been evaluated in various analyses over the last decade. In Ref. [50] the leading Yukawa 
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corrections to A/H —)■ X 2 X 2 ; xtXi the rMSSM have been evalnated, employing an on-shell 
(OS) scheme (referring to Ref. [51], bnt withont providing further details). Next, in Ref. [52] 
the full one-loop corrections to A/H ^ xtXc' — 1)2) have been presented in the 

rMSSM (again without details about the OS like scheme). An effective Lagrangian approach 
for heavy neutral Higgs boson decays in the rMSSM was published in Ref. [53]. The full 
one-loop corrections to all heavy Higgs decays to charginos and neutralinos in the rMSSM 
in the DR scheme was published in the code HFOLD M. More recently also evaluations 
of Higgs boson decays to charginos and neutralinos in the cMSSM became available. In 
Ref. [51] the decays hi —)> xtx^o = 1)2,3; c, c' = 1,2) were presented, together with a 
short discussion of different renormalization schemes (see Sect. 12 .21) and brief analysis of the 
dependence on the phases of /i. Mi and the trilinear Higgs stop coupling, At. The decays 
hi — )■ Xn Xn' (* = 1)2,3; n,n' = 1,2, 3,4) were calculated in Ref. [55], where the numerical 
analysis concentrated on h 2,3 ^ X 2 X 2 dependence on the phase of At. The latter 

two references are close to the calculations presented in this paper. Small differences in the 
renormalization in the chargino/neutralino sector exist (see Sect. 12.21 and Ref. [S]), where 
we use consistently the scheme detailed in Ref. [15] for all two-body decays simultaneously. 
In our numerical analysis we focus on the one parameter with a possible complex phase 
entering at the tree-level. Mi, see the discussion in Sect. El A short numerical comparison 
with the literature, in particular with Refs. [10l[511[55], will be given in Sect. 14.11 

In this paper we present a full one-loop calculation for all two-body decay channels of 
the Higgs bosons into charginos and neutralinos in the cMSSM, taking into account soft 
and hard QED radiation. In Sect. [2] we briefly review the relevant sectors of the cMSSM. 
Details about the calculation can be found in Sect. E] and the numerical results for all decay 
channels are presented in Sect. |4] (including comments on comparisons with results from other 
groups). The conclusions can be found in Sect. El It is planned to implement the evaluation 
of the branching ratios of the heavy Higgs bosons into the Fortran code FeynHiggs [26II30] . 
together with an automated choice of the renormalization scheme valid for the full cMSSM 
parameter space. 

2 The complex MSSM 

The channels CO - o are calculated at the one-loop level, including soft and hard QED 
radiation. This requires the simultaneous renormalization of several sectors of the cMSSM: 
the Higgs and gauge boson sector as well as the chargino/neutralino sector. In the following 
subsections we very briefly review these sectors and their renormalization. 

2.1 The Higgs- and gauge-boson sector 

The Higgs- and gauge-boson sector follow strictly Ref. [45] and references therein (see es¬ 
pecially Ref. 1251). This defines in particular the counterterm 5 tan {3 = 6ty, as well as the 
counterterms for the Z boson mass, (5M|, and for the sine of the weak mixing angle, 5sw 
(with Sw = -^/l — — M^/M|, where Mw denotes the W boson mass). 
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2.2 The chargino/neutralino sector 

The chargino/neutralino sector is also described in detail in Ref. |15] and references therein, 
see in particular Refs. [161 - H8] . In this paper we use the so called “CCN” scheme, i.e. OS 
conditions for two charginos and one neutralino, which we chose to be the lightest one. 
Renormalizing the two charged states OS, i.e. ensuring that they have the same mass at the 
tree- and at the loop-level is (in general) crucial for the cancellation of the IR divergencies. 
In the notation of Ref. [IS] we used: 

SinoScheme = CCN[1] hxed CCN scheme with on-shell Xi • 

This dehnes in particular the counterterm where denotes the Higgs mixing parameter. 
This scheme yields numerically stable results for \Mi\ < IM 2 I, \fi\ and M 2 7 ^ /i, i.e. the lightest 
neutralino is bino-like and dehnes the counterterm for Mi [15IH5] . In the numerical analysis 
this mass pattern holds. Switching to a different mass pattern, e.g. with IM 2 I < \Mi \ and/or 
M 2 ~ /i requires to switch to a different renormalization scheme [^5E^ . While these schemes 
are implemented into the FeynArts/FormCalc framework [l5|, so far no automated choice of 
the renormalization scheme has been devised. For simplicity we stick to the CCN[1] scheme 
with a matching choice of SUSY parameters, see Sect. 14.21 

Since both chargino masses m-± and the lightest neutralino mass m^o have been chosen 

as independent parameters, the one-loop masses of the heavier neutralinos Xn = 2,3,4) 
are obtained from the tree-level ones via the shifts [55] 



where the renormalization constants 6Z and 6Z can be found in Ref. [45]. For all externally 
appearing neutralino masses we use the (shifted) “on-shell” masses: 

= m^o + Am^o . (5) 

Xn 2Cn An v / 

In order to yield UV-hnite results we use the tree-level values m^o for all internally appearing 
neutralino masses in loop calculations. 

2.3 The fermion/sfermion sector 

To be in accordance with Ref. [42], we use shifted (s)fermion masses in the loop corrections. 
As requirement for these shifts one needs the renormalization of the fermion/sfermion sector: 

• The renormalization of the fermion sector is described in detail in Ref. [45] and refer¬ 
ences therein. For simplihcation we use the DR renormalization for all three generations 
of down-type quarks and leptons, again in the notation of Ref. [45] : 

UVMfl[4, _] = UVDivergentPart DR renormalization for rud, mg, 

UVMfl[2, _] = UVDivergentPart DR renormalization for rUe, m^, 
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• The renormalization of the sfermion sector differs slightly from the one described in 
Ref. [IS]. For the squark sector we follow Refs. |131lll| (which agrees with the renor¬ 
malization scheme used in Refs. [ISIIIZIESIISZI) and the slepton sector can be found 
in Ref. [12]. Concerning our notation we denote as ^Qg igUg bg Eg^ “diagonal” 
soft SUSY-breaking parameters for the SU{2) squark, slepton doublet, the u-, d-type 
squark singlet, and the e-type slepton singlet, respectively, where g is the generation 
index. Furthermore we use d g for tho trilinear Higgs-scalar u-, d-, e-type fermion 
couplings, respectively. 


3 Calculation of loop diagrams 

In this section we give some details about the calculation of the higher-order corrections 
to the partial decay widths of Higgs bosons. Sample diagrams for the decays hi —)■ XcXt' 
{i = 1,2,3; c,c' = 1,2), hi xlxl' (* = 1,2,3; n,n' = 1,2, 3,4) and H+ xlxt = 
1, 2 , 3,4; c = 1, 2 ) are shown in Figs. [Tl[2]and[3l respectively. Not shown are the diagrams for 
real (hard and soft) photon radiation. They are obtained from the corresponding tree-level 
diagrams by attaching a photon to the electrically charged particles. The internal generically 
depicted particles in Figs. [H [2] and [S] are labeled as follows: F can be a SM fermion /, 
chargino neutralino S can be a sfermion fs or a Higgs (Goldstone) boson hi (G); V 
can be a photon 7 or a massive SM gauge boson, Z or W^. For internally appearing Higgs 
bosons no higher-order corrections to their masses or couplings are taken into account; these 
corrections would correspond to effects beyond one-loop order|i| For external Higgs bosons, 
as discussed in Ref. [26], the appropriate Z factors are applied and OS masses (including 
higher-order corrections) are used [26], obtained with FeynHiggs [261130] . 

Also not shown are the diagrams with a Higgs boson-gauge/Goldstone self-energy con¬ 
tribution on the external Higgs boson leg. They appear in the decay hi -b xFX)!', Fig- [Hand 
hi —)■ XnXn'y Fig. 121 with a hi-ZjG transition and in the decay —)> Fig. [31 with a 

//±_VP±/( 7 ± transition! 

Furthermore, in general, in Figs. [H^E] we have omitted diagrams with self-energy type 
corrections of external (on-shell) particles. While the contributions from the real parts of 
the loop functions are taken into account via the renormalization constants defined by OS 
renormalization conditions, the contributions coming from the imaginary part of the loop 
functions can result in an additional (real) correction if multiplied by complex parameters. 
In the analytical and numerical evaluation, these diagrams have been taken into account via 
the prescription described in Ref. [45] . 

Within our one-loop calculation we neglect finite width effects that can help to cure 
threshold singularities. Gonsequently, in the close vicinity of those thresholds our calculation 
does not give a reliable result. Switching to a complex mass scheme [3S] would be another 
possibility to cure this problem, but its application is beyond the scope of our paper. 

The diagrams and corresponding amplitudes have been obtained with FeynArts [59] . 
The model hie, including the MSSM counterterms, is largely based on Ref. [45], however 

^ We found that using loop corrected Higgs boson masses in the loops leads to a UV divergent result. 

^ From a technical point of view, the H^-W ^/transitions have been absorbed into the respective 
counterterms, while the hi-ZjG transitions have been calculated explicitly. 
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Figure 1: Generic Feynman diagrams for the decay hi —)■ XcXt' (* = 1 ; 2 , 3; c, c' = 1, 2). F can 
be a SM fermion, chargino, neutralino; S can be a sfermion or a Higgs/Goldstone boson; V can 
be a 7 , Z, Not shown are the diagrams with a h—Z or h—G transition contribution on the 
external Higgs boson leg. 



Figure 2: Generic Feynman diagrams for the decay hi —)■ XnXn' (* ~ 15 2 , 3 ; n,n' = 1 , 2 , 3 , 4 ). 
F can be a SM fermion, chargino, neutralino; S can be a sfermion or a Higgs/Goldstone boson; 
V can be a Z or W^. Not shown are the diagrams with a hi-Z or h—G transition contribution 
on the external Higgs boson leg. 


adjusted to match exactly the renormalization prescription described in Sect. El The further 
evaluation has been performed with FormCalc and LoopTools [60] . 

Ultraviolet divergences 

As regularization scheme for the UV divergences we have used constrained differential renor¬ 
malization [61], which has been shown to be equivalent to dimensional reduction [62] at the 
one-loop level [60]. Thus the employed regularization scheme preserves SUSY [HSlEl] and 
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Figure 3: Generic Feynman diagrams for the decay XnXt = 1)2, 3,4; c = 1,2). (It 

should be noted that all arrows are inverted in case of a H~ decay.) F can be a SM fermion, 
chargino, neutralino; S can be a sfermion or a Higgs/Goldstone boson; V can be a 7, Z, W^. 
Not shown are the diagrams with a or H^-G^ transition contribution on the external 

Higgs boson leg. 


guarantees that the SUSY relations are kept intact, e.g. that the gauge couplings of the 
SM vertices and the Yukawa couplings of the corresponding SUSY vertices also coincide to 
one-loop order in the SUSY limit. Therefore no additional shifts, which might occur when 
using a different regularization scheme, arise. All UV divergences cancel in the hnal result. 

Infrared divergences 

The IR divergences from diagrams with an internal photon have to cancel with the ones from 
the corresponding real soft radiation. In the case of QED we have included the soft photon 
contribution following the description given in Ref. |ST]. The IR divergences arising from 
the diagrams involving a 7 are regularized by introducing a photon mass parameter, A. All 
IR divergences, i.e. all divergences in the limit A —?■ 0, cancel once virtual and real diagrams 
for one decay channel are added. 

Tree-level formulas 

For completeness we show here also the formulas that have been used to calculate the tree- 
level decay widths: 


"'tree/ 


XnXf) = X^)lP + \C{H^, - m|o - m|±) 


- 4 Re Xn, &L xl, rn^i 

-^ c= 1,2; n = 1,2,3,4 , 

16 7rM^± ^ ’ ’ ’ ’ ’ 


X 


( 6 ) 
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ntree 


{hi Xc xt') = 


{\c{hi, x^ , xj)Lr+ \c{hi, Xe, xt')Rn - ^1+) 

- 4 Re {C{hi, x~, xt')! C{hi, x~, m^+ 

(z = l,2,3;c,c' = l,2), 


16 vrm^. 


-'tree/ 


{hi^xlxl')= {\C{hi,Xn,Xn')L\^+ \C{hi,Xn,Xn>)Rn {T^hi-rnlo -mlo ] 

A,„ ^TJ' 


4 Re {C{hi, Xn , xI'Yl C{hh xl , XnOij} 


X 


A 1 / 2(^2 2 ^^2 ) 


16 vrm^^ 


(i = 1, 2, 3; n, n' = 1, 2, 3,4), 


( 7 ) 


( 8 ) 


where X{x,y,z) = {x — y — zY — iyz. The couplings C{a,b,c) can be found in the 
FeynArts model hies, see Ref. [6S]. C{a,b,c)L^R denote the part of the coupling which 
is proportional to (1 =F 75)/2. For the later interpretation of the results in the numerical 
analysis the following should be kept in mind. In Eqs. ([7]), (|H]) the couplings of the Higgs to 
charginos/neutralinos result in a relative plus (minus) sign between the two terms (in the 
hrst and second line of each equation, respectively) for hi being a CP-odd (CP-even) Higgs, 
leading to an enhancement (suppression) of the decay width. In case of equal hnal state 
masses (c = c' or n = n') one hnds in FeynArts convent iorU C(hj, x, x)l = ~C{hi,x,x)*j^ =: 
C{hi, X, x), and the general structure of the tree-level decay width simplihes for real param¬ 
eters to 


^ xx) 
r^^^Yh, ^ xx) 


\C{hi,x,x)\^ 

8 TT 

\C{hi,X,x)\^ 

SttuiY 


hi 



for hi CP-odd, 

(9) 


for hi CP-even. 

(10) 


The latter decay width exhibits a p-wave suppression. 


4 Numerical analysis 

In this section we present the comparisons with results from other groups and our numerical 
analysis of the light and heavy Higgs boson decay channels into charginos and neutralinos in 
the cMSSM. In the various hgures below we show the partial decay widths and their relative 
correction at the tree-level (“tree”) and at the one-loop level (“full”). 

4.1 Comparisons 

We performed exhaustive comparisons with results from other groups for Higgs boson decays 
into charginos and neutralinos. Most of these comparisons were restricted to the MSSM with 
real parameters. 

^ It should be noted that the convention for Feynman rules in Ref. [2] and in Ref. [65] differ by a global 
factor of —i, which would formally lead to C{hi,x,x)L = +C{hi,x,x)R- However, the physics outcome 
remains, of course, unchanged. 











Table 1: Comparison of the 1-loop corrected partial decay widths (in GeV) with Ref. (5^ . 


Process 

Mass 

Ref. [52] 

Ref. [55] 

FeynTools 

71° ^ xtx/ 

0 

0 

0.85 

0.80 

0.83 

71° ^ xtxi 

Ma = 800 

1.00 

0.91 

0.96 


Mh = 800 

0.63 

0.58 

0.64 

^° ^ x^xr 

Mh = 900 

0.73 

0.70 

0.75 


• A comparison with Ref. [50] (in the rMSSM) gave an overall qualitative agreement for 
the decays H/A xtXiy was to be expected, because that work took into account 
only the leading Yukawa corrections and used a different renormalization scheme. On 
the other hand, we omit a comparison with the results for H/A ^ X2X2 of Ref. |50] . 
because in their set-up neutralino masses were used as input parameters, which is 
rather difficult to adapt to our numerical analysis set-up. 

• For the comparison with Ref. [52] we calculated the decays A/H ^ xtXi 0{a{Mz)), 
using their input parameters as far as possible. We found good (qualitative) agreement 
with Ref. [52] (where the calculation was restricted to the rMSSM). We successfully 
reproduced their Figs. 2, 4, 5 and 6, where only a small difference remains due to the 
different renormalization schemes, see also Tab. [Tjwith differences below 5%. 

• We performed a numerical comparison with Ref. [53] for the decay H/A —)■ xfx/ at 
0{a{Mz)) (in the rMSSM), see Tab. d] where the columns for Ref. [52] and Ref. [55] 
where taken over from Tab. 1 of the hrst article in Ref. [55]. Their set-up differs from 
ours in the renormalization of the chargino/neutralino sector, leading to different loop 
corrections. Furthermore they used an “effective one-loop Lagrangian”. Nevertheless, 
using their input parameters as far as possible, we found differences below the 10% level. 

• A numerical comparison with the program HFOLD Ref. [ID] at the benchmark point 
SPSla‘ (proposed in the SPA project [66]) can be found in Tab. [21 Only for this point 
sufficient details about the chargino/neutralino masses was available for a numerical 
comparison. In Tab. [2] we show the full one-loop results of HFOLD, using DR masses 
for the internal and external particles, corresponding to the full DR renormalization 
used in the code (where the renormalization scale was set to 1 TeV [66]). Our results, 
labeled FeynTools, are evaluated using our renormalization scheme, but inserting the 
HFOLD DR masses. In the tree-level results we hnd more than 10 digits agreement and 
in the full results we hnd agreement of 3% - 15% (7% on average). 

HFOLD also offers to switch to (the recommended) OS masses for the external particles. 
In this case, we are including in our calculation m7o as described in Sect. 12.21 Eq. ([5]), 
but using the same OS Higgs boson masses as in HFOLD. With it the agreement between 
the two calculations is 11% on average (5% - 22%). 

• Decays of ^2,3 to charginos in the cMSSM at the full one-loop level have been nu- 
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Table 2: Comparison of the 1-loop corrected partial decay widths (in 10 ^ GeV) with HFOLD. 


Process 

OS masses 

DR masses 

HFOLD FeynTools 

HFOLD FeynTools 


0.1381 

0.1648 

0.1046 

0.1229 


0.4584 

0.4908 

0.2690 

0.2828 


0.2061 

0.2259 

0.0117 

0.0111 


0.5262 

0.5672 

0.0345 

0.0332 

71° ^ X?X? 

0.2044 

0.2404 

0.1704 

0.2016 

71° ^ X?X° 

0.9693 

1.0248 

0.7334 

0.7750 

71° ^ X°X° 

1.1652 

1.0747 

0.3966 

0.3791 

71° ^ xtxl 

2.8604 

2.6454 

1.0236 

0.9928 

^ X?X^ 

1.2981 

1.4307 

0.9333 

0.9996 

^ X^X^ 

0.0063 

0.0081 

0.0026 

0.0030 


merically compared with Ref. [5l] using their latest FeynArts/FormCalc model file 
implementation. We found overall agreement better than 1% in the loop corrections 
for real and complex parameters]^ 


• /i2,3 boson decays into X 2 X 2 ^^e cMSSM have been analyzed in Ref. [55]. Again we 
had to use here the latest FeynArts/FormCalc model hie implementation of Ref. [5l| 
(which bases mainly on code from Ref. [55]) for the same reasons as described in the 
previous item. In comparison with that model hie [M] we found overall agreement 
better than 2% in the loop corrections for real and complex parameters. 


4.2 Parameter settings 

The renormalization scale fiR has been set to the mass of the decaying Higgs boson. The 
SM parameters are chosen as follows; see also [5B] : 

• Fermion masses (on-shell masses, if not indicated diherently): 


rUe = 0.510998928 MeV , 
= 105.65837515 MeV, 
nir = 1776.82 MeV, 
rriu = 68.7 MeV , 
rric = 1.275 GeV, 
mt = 173.21 GeV, 


= 0 MeV , 
rriui, = 0 MeV , 

= 0 MeV , 
rrid = 68.7 MeV , 
rris = 95.0 MeV , 
mfe = 4.18 GeV . 


( 11 ) 


^ It should be noted that the original code used for Ref. [54] is no longer available [62], where we found 
some numerical differences with the results shown in Ref. [^ in the case of complex parameters. 


10 









According to Ref. [68], rris is an estimate of a so-called ’’current quark mass” in the 
MS scheme at the scale /i ~ 2 GeV. rric = mc{mc) and = mb{mf,) are the ’’running” 
masses in the MS schemeJl and are effective parameters, calculated through 
the hadronic contributions to 

= ^ E ~0'027723. (12) 

f=u,c,d,s,b f 

• Gauge boson masses: 

Mz = 91.1876 GeV, Mw = 80.385 GeV. (13) 


• Goupling constant: 


a(0) = 1/137.0359895. (14) 

The Higgs sector quantities (masses, mixings, etc.) have been evaluated using FeynHiggs 
(version 2.10.2) [261130] . 

We emphasize again that the analytical calculation has been done for all decays into 
charginos/neutralinos. Results are shown for some representative numerical examples. The 
parameters are chosen according to the scenarios, SX (X = 1,2,...,5), shown in Tab.[3l unless 
otherwise noted. The scenarios are dehned such that a maximum number of decay modes are 
open simultaneously to permit an analysis of all channels, i.e. not picking specihc parameters 
for each decay. For the same reason we do not demand that the lightest Higgs boson has a 
mass around ~ 125 GeV, although for most of the parameter space this is given. For the 
light Higgs we will show the variation with Mh±, |/i|. Mi and (where the latter denotes 
the phase of the gaugino mass parameter Mi), whereas for the heavy Higgs bosons we will 
analyze the variation of Mh± and (pMi- 

The numerical results shown in the next subsections are of course dependent on choice of 
the SUSY parameters. Nevertheless, they give an idea of the relevance of the full one-loop 
corrections. Ghannels (and their respective one-loop corrections) that may look unobservable 
due to the smallness of their decay width in the plots shown below, could become important 
if other channels are kinematically forbidden. 

4.3 Full one-loop results for varying Mi, and (fMi 

The results shown in this and the following subsections consist of “tree”, which denotes 
the tree-level value and of “full”, which is the partial decay width including all one-loop 
corrections as described in Sect. [3l We restrict ourselves to the analysis of the decay widths 
themselves, since the one-loop effects on the branching ratios are strongly parameter depen¬ 
dent, as discussed in the previous subsection. 

When performing an analysis involving complex parameters it should be noted that the 
results for physical observables are affected only by certain combinations of the complex 
phases of the parameters /i, the trilinear couplings Af and the gaugino mass parameters 

® It should be noted, that in the analysis below, we use the DR mass from Eq.(19) of Ref. [42] . 
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Table 3: MSSM parameters for the initial numerical investigation; all parameters (except of 
tjs) are in GeV (calculated masses are rounded to 1 MeV). In our analysis 
Mi^ and are chosen such that the values of mi^, mi^, rriy^ and are realized. For 
the sfermion sector the shifts in Mq jj{dg) and Mi^{eg) as defined in Ref. [42] are taken into 
account. The values for the trilinear sfermion Higgs couplings, {g = 1,2,3; identical for 

all g) are chosen such that charge- and/or color-breaking minima are avoided [65]- It should be 
noted that for the first and second generation of sfermions we chose instead ^ = 1500 GeV 
and Mgfj = 2000 GeV. For the neutralino sector the shifts in Eq. (HI) are taken into account. 


Seen. 

tp 

h 


Ad 

^9 

Ae 

|Mi| 

Ms 

Mg 


mia ^b2 rnf2 

SX 

10 

500 

1200 

600 

1000 

300 

600 

1500 

394 

771 582 280 309 


Seen. 

Mh± 

SI 

700 

S2 

900 

S3 

1000 

S4 

1200 

S5 

1400 

rrihi 

123.487 

123.509 

123.517 

123.529 

123.539 

mh2 

694.483 

895.594 

996.769 

1195.095 

1397.300 

rriha 

695.425 

896.931 

996.818 

1197.407 

1398.600 

rrif^ 

282.705 

282.573 

282.517 

282.420 

282.336 


513.733 

513.621 

513.578 

513.509 

513.455 

m^o 

295.269 

295.269 

295.269 

295.269 

295.269 

m^o 

A2 

476.772 

476.763 

476.759 

476.753 

476.748 

m^o 

A3 

496.992 

496.988 

496.986 

496.983 

496.980 


632.326 

632.324 

632.324 

632.323 

632.322 


472.534 

472.534 

472.534 

472.534 

472.534 

^xt 

632.167 

632.167 

632.167 

632.167 

632.167 


4 ^1, 2,3 [ZOllTT]. It is possible, for instance, to rotate the phase ipM 2 away. Experimental 
constraints on the (combinations of) complex phases arise, in particular, from their contri¬ 
butions to electric dipole moments of the electron and the neutron (see Refs. [721173] and 
references therein), of the deuteron [71| and of heavy quarks [75]. While SM contributions 
enter only at the three-loop level, due to its complex phases the MSSM can contribute al¬ 
ready at one-loop order. Large phases in the hrst two generations of sfermions can only be 
accommodated if these generations are assumed to be very heavy ug or large cancellations 
occur see, however, the discussion in Ref. [78]. A review can be found in Ref. [79] . 
Accordingly (using the convention that ipMa = 0, as done in this paper), in particular, the 
phase (pg is tightly constrained [80] , while the bounds on the phases of the third generation 
trilinear couplings are much weaker. Setting ipg = 0 and pAf = 0 leaves us with as the 
only complex valued parameter. 

Since now the complex gaugino mass parameter Mi can appear in the couplings, con¬ 
tributions from absorptive parts of self-energy type corrections on external legs can arise. 
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The corresponding formulas for an inclusion of these absorptive contributions via hnite wave 
function correction factors can be found in Refs. [HIHS] , 

We begin the numerical analysis with partial decay widths of evaluated as a function 
of Mh±, starting at Mh± = 600 GeV up to Mh± = 1-6 TeV, which roughly coincides with 
the reach of the LHC for high-luminosity running as well as an e’''e“ collider with a center- 
of-mass energy up to ^/s ~ 3 TeV [16]. Then we turn to the hi {i = 1, 2, 3) decays. Finally, 
it should be noted that we expect from the tree-level formulas Eqs. ([S]) ^ dH]) that the decay 
widths increase (roughly) linearly with the corresponding Higgs boson masses. 

4.3.1 decays into charginos/neutralinos 

In Figs, m - [11] we show the results for the processes —)■ ~ 1)2, 3,4; c = 1,2) 

as a function of Mu± and as a function of the relevant complex phase ^pmi- These are of 
particular interest for LHC analyses [8T1[82] (as emphasized in Sect. [1]). The various visible 
(or hardly visible) dips/thresholds occurring for different values of Mu± in the plots are 
summarized in Tab. 01 labeled TCI to TC7. 

We start with the decay Xixf- fh® plof of Fig. 0] the hrst (small) dip is 

the threshold TCI, see Tab. 0] The second (large) dip is an effect due to the threshold TC2. 
The third “apparently single” dip is in reality two dips coming from the thresholds TC3 and 
TC4. The fourth (small) dip is the threshold TC5 and the last (large) one is the threshold 
TC6. The size of the corrections of the partial decay widths can be especially large very 
close to the production thresholci^ from which on the considered decay mode is kinematically 
possible. Away from this production threshold relative corrections of ~ -|-10% are found. 

In the right plot of Fig. 0] we show the results for the complex phase (fMi varied for 
Mh± = 1000 GeV. The full corrections are up to ~ -|-13% at ipMi = 180°. At (fMi = 90° the 
{H~) full corrections reach ~ -|-12% (~ 4-10%). Clearly visible is the CP-asymmetry 
for the decays of the and H~, which can reach the level of several per-cent. 

In Fig. [5]we show the results for —)■ Xixf- The tree-level decay width F(iL^ —)■ Xixt) 

is accidently very small for the parameter set chosen, see Tab. [S] Because of this small¬ 
ness, the relative size of the one-loop correction becomes larger than the tree-level result, 
and can even turn the decay width, oc |A4tree|^ + 2RejAltj-eenegative. There¬ 
fore, in this case we added |A4i_ioopP to the full one-loop result to obtain a positive de¬ 
cay width. In the left plot the hrst (large) spike is the threshold TC2, see Tab. 0] en¬ 
hanced through the two-loop contribution |A4i_ioopP (i-e. without the explicit two-loop 
correction the spike would be a “usual dip”). The second “apparently single spike” 
(hardly visible) is (again) in reality the two thresholds TC3 and TC4. The next (ap¬ 
parently single) “dip” is in reality two steps (anomalous thresholds, see Ref. [SB]) traced 
back to the C-functions ) at Mh± ps 1126 GeV and 

’ ’ ^ X2 Xi X4 Xi |— 

Co,i, 2 (Tf|r±, m?±, m- 0 ,m- 0 , Af^) at Mh± ~ 1129 GeVhl Not visible (in the plot) is a 

® It should be noted that a calculation very close to the production threshold requires the inclusion of 
additional (nonrelativistic) contributions, which is beyond the scope of this paper. Consequently, very close 
to the production threshold our calculation (at the tree- and loop-level) does not provide a very accurate 
description of the decay width. 

^ In addition both steps are contorted through the higher order contributions |AIi-ioopp. 
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spike, which is the threshold TC5. The last spike is the threshold TC6. Relative corrections 
of ~ —90% are found at Mh± = 1000 GeV (see Tab. [3]), where it should be kept in mind 
that the tree-level is already accidentally very small and thus loop corrections can have a 
relatively large impact. 

In the right plot of Fig. [5]the results are shown for S3 as a function of At ipMi = 180° 
the full corrections reach ~ —55%, again related to the accidentally small tree-level result. 
At ipMi = 90° the {H~) full corrections reach ~ —59% (~ —62%), showing a small 
CP-asymmetry. 

Next, in Fig. Elthe results for X 2 xt displayed. In the left plot the results are 

shown as a function of Mh±. The four visible dips here are exactly the same as in Fig. |4] 
(described above), beginning at Mh± = 976 GeV. Relative corrections of ~ -|-33% (~ -|-21%) 
are found at Mh± = 1000 GeV {Mh± = 1400 GeV), see Tab. [3l 

In the right plot the results are displayed as a function of (fMi in S3, i.e. for Mh± = 
1000 GeV. One can see that the size of the tree-level as well as the corrections to the partial 
decay width vary substantially with the complex phase tpMi ■ For all the full corrections 
lie between -|-29% and -|-70%G At (fMi = 90° the {H~) full one-loop corrections reach 
~ -|-29% (~ 4-48%), i.e. the CP-asymmetry is rather large with ~ ±19%. 

The decay is shown in Fig. [3 In the left plot the results are shown as a 

function of Mu±. The hrst (hardly visible) dip is (again) the threshold TG5, see Tab. [Hand 
the second (large) one is the threshold TG6. The decay width turns out to be relatively 
large at 0{1 GeV). Relative corrections of ~ ±6% (~ ±4%) are found at Mh± = 1200 GeV 
{Mh± = 1400 GeV), see Tab. |3l 

In the right plot of Fig. [3 the results are displayed as a function of (pMi in S4, i.e. for 
Mh± = 1200 GeV. The full corrections at ifMi = 180° reach ~ ±7%. On the other hand 
it can be seen that the variation with and the CP-asymmetries (~ ±0.1%) are rather 

small. 

Next, in Fig. [8] the results for —)> XsXi' nre displayed. In the left plot the results are 
shown as a function of Mh±. Here the four visible dips are the same as in Fig. 01 beginning at 
Mh± = 976 GeV. Relative corrections of ~ —18% (~ —10%) are found at Mh± = 1000 GeV 
= 1400 GeV), see Tab. [3 

In the right plot the results are displayed as a function of ipMi in S3. One can see that 
again the tree-level results as well as the size of the corrections to the partial decay width 
vary substantially with the complex phase ifMi ■ The full corrections can reach ~ —24% and 
the CP-asymmetry is found to be small at the level of ~ ±2%. 

In Fig. Owe show the results for xtxf- In the left plot (with Mh± varied) the dip 

is (again) the threshold TG6, see Tab. 01 The decay width is found to be of the same order as 
for —)■ xlxt- One-loop corrections of ~ ±5% (~ ±4%) are found at Mh± = 1200 GeV 

in S4 {Mh± = 1400 GeV in S5), see Tab. O 

In the right plot of Fig. Othe results are shown for S4 as a function of At ipMi = 180° 
the full corrections reach ~ ±5%. At = 90° the {H~) full corrections reach ~ ±6% 
(~ ±5%), i.e. the CP-asymmetries are at the level of ~ ±1%. 

® It should be noted that the loop corrections can reach -1-70% of the tree results because at ^Mi = 180° 
the tree-level decay width is accidently small, see the right plot of Fig. 01 
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Table 4: Thresholds in charged Higgs boson decays. 


TGI 

Mh± ~ 907 GeV 

TG2 

Mh± = 976 GeV 

TG3 

Mh± ~ 1105 GeV 

TG4 

Mh± ~ 1108 GeV 

TG5 

Mh± ~ 1135 GeV 

TG6 

Mh± ~ 1284 GeV 

TG7 

Mh± ~ 1353 GeV 


^k + 

m~± + m^o = Mh± 

Xi X4 

m^± + m^o = Mh± 

^xt + ^xl = 

'^k + 

'^k + 


We finish the charged Higgs-boson analysis with the decays involving the heaviest neu- 
tralino in Figs. ITO] and fTTl showing the results for r(ih=*= ^ XaX^) and T{H^ —>■ 
respectively. 

In the left plot of Fig. [10] the hrst dip (not visible in the plot) is the threshold TC4, see 
Tab. m The second (small) dip is (again) the threshold TC5 and the third (large) dip is the 
threshold TC6. The hrst step (not visible in the plot) at Mu± k, 1136 GeV is the anomalous 
threshold of the G-functions Go,i,2(Tf^±, m^±, m-o, m-o,'^-±, Tf|). The second anomalous 

threshold at Mh± ~ 1340 GeV is caused by Go,i, 2 (^^|o, m|). The last 

dip (also not visible) is the threshold TG7. The decay width is again found at 0{1 GeV) 
with relative corrections of ~ +6% in S4 (see Tab. [3l). 

In the right plot of Fig. [10] we show the complex phase (p^i varied at Mh± = 1200 GeV. 
The full corrections are up to ~ +6% at = 180°. Here the asymmetries are extremely 
small and hardly visible. 

Finally, we discuss the decay X 4 X 2 Fig. [TT] The overall size of this decay width 

(with real phases) is (accidentally) very small around 1 x 10“^ GeV. Gonsequently, the loop 
corrections, can be larger than the tree-level result. In the left plot the results are shown as 
a function of Mh±- The (small) dip is the threshold TG6, see Tab. ]4] Relative corrections 
of ~ -|-56% are found at Mh± = 1400 GeV (see Tab. [3]). 

In the right plot the results are displayed as a function of tpMi in S5. One can see that the 
size of the corrections to the partial decay width vary substantially with the complex phase 
9 ?Ml- For all 9 ?Ml the full corrections deviate between -|-40% and -|-146%. (The latter value 
is reached aX ip= 180° where the tree is extremely small ~ 1 x 10-^ GeV.) At <^Mi = 90° 
the {H~) full one-loop corrections reach ~ -|-40% (~ -|-103%), i.e. the CP-asymmetries 
are very large with ~ ±60%. 

Overall, for the charged Higgs boson decays to a chargino/neutralino pair we observe, as 
expected, an increasing decay width oc Mm±|j The full one-loop corrections reach a level of 
10% for decay widths of (9(1 GeV). The variation with pMi is found largest for very small 
decay widths, but can reach the level of 10 — 50% for widths at or below the 1 GeV level. 
The CP-asymmetries exceed the level of a few per-cent only for very small decay widths. 


An exception are the loop corrections in the left plot of Fig. [SI because there we added |Afi-ioop| 
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r/GeV r/GeV ^ ^0;^± 



Mh± 


‘PMi 


Figure 4: T{H^ —)■ Xixf)- Tree-level and full one-loop corrected partial decay widths are 
shown. The left plot shows the partial decay width with Mh± varied. The right plot shows the 
complex phase ipMi varied with parameters chosen according to S3 (see Tab. [3]). 


rxlO-^GeV rxlO-VGeV ^ 



Mh± 


V^Mi 


Figure 5: T{H^ —)■ Xixf). Tree-level and full one-loop corrected partial decay widths are 
shown. The left plot shows the partial decay width with Mh± varied. The right plot shows the 
complex phase varied with parameters chosen according to S3 (see Tab. [3]). 
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r/GeV £f± ^ r/GeV ^ ;^o^± 



Figure 6: T{H^ —)■ xlxf)- Tree-level and full one-loop corrected partial decay widths are 
shown. The left plot shows the partial decay width with Mh± varied. The right plot shows the 
complex phase ipMi varied with parameters chosen according to S3 (see Tab. [3]). 


r/GeV ^ x°2X2 r/GeV ff± ^ 





Figure 7: T{H^ —)■ X 2 X 2 )■ Tree-level and full one-loop corrected partial decay widths are 
shown. The left plot shows the partial decay width with Mh± varied. The right plot shows the 
complex phase varied with parameters chosen according to S4 (see Tab. [3]). 
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r/GeV i/± ^ ;^o^± r/GeV ^ 



Figure 8: T{H^ —)■ Xaxf)- Tree-level and full one-loop corrected partial decay widths are 
shown. The left plot shows the partial decay width with Mh± varied. The right plot shows the 
complex phase ipMi varied with parameters chosen according to S3 (see Tab. [3]). 


r/GeV ^ r/GeV ff± ^ 







Figure 9: T{H^ —)■ X 3 X 2 )■ Tree-level and full one-loop corrected partial decay widths are 
shown. The left plot shows the partial decay width with Mh± varied. The right plot shows the 
complex phase varied with parameters chosen according to S4 (see Tab. [3]). 
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r/GeV ^ r/GeV 



Mh± 


‘PMi 


Figure 10: T{H^ X^xf)- Tree-level and full one-loop corrected partial decay widths are 
shown. The left plot shows the partial decay width with Mh± varied. The right plot shows the 
complex phase ipMi varied with parameters chosen according to S4 (see Tab. [3]). 


rxlO-^GeV rxlO-VGeV xlX^ 



Figure 11: T{H^ —)■ XiXt)- Tree-level and full one-loop corrected partial decay widths are 
shown. The left plot shows the partial decay width with Mh± varied. The right plot shows the 
complex phase 99^1 varied with parameters chosen according to S5 (see Tab. [3]). 
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4.3.2 hi decays into charginos and neutralinos 

We now turn to the decay modes hi XcXt' (* = 2,3; c, c' = 1,2) and hi —)■ 

{i = 1, 2, 3; n, n' = 1, 2, 3,4). Results are shown in the Figs. [T2]-[25l 

Before discussing every hgure in detail, it should be noted that there is a subtlety con¬ 
cerning the mixture of the hi bosons. Depending on the input parameters, the higher-order 
corrections to the three neutral Higgs boson masses can vary substantially. The mass order¬ 
ing < mh 2 < (as performed automatically by FeynHiggs), even in the case of real 
parameters, can yield a heavy CP-even Higgs mass higher or lower than the (heavy) CP-odd 
Higgs mass. Such a transition in the mass ordering (or “mass crossing”) is accompanied by an 
abrupt change in the Higgs mixing matrix Z0 For our input parameters (see Tab. [3]) there 
are two (possible) crossings. The hrst (called “MCI” below) appears at ~ 1006 GeV. 
Before the crossing we hnd h^ ^ H {h^ ~ A), whereas after the crossing it changes to 
h 2 ~ H (hs ~ H). The second crossing (called “MC2”) is found at Mh± ~ 1532 GeV, i.e. 
the changing of the mixture from ^2 ~ H {h^ ~ H) to h 2 ~ P {h^ ~ A). Very close to 
the mass crossings the Z matrix can yield small numerical instabilities. As an example, for 
1532 GeV < Mh± < 1536 GeV the Z matrix causes structures appearing similar to “usual” 
dips from thresholds (see also the discussion in Ref. |12]). All the dips/thresholds (some are 
hardly visible) appearing in the hgures below are listed in Tab. [5], labeled as TNI to TN13. 

hi decays into charginos 

In this subsection we analyze the decays of the heavy neutral Higgs bosons into charginos. 
We start with the decay hi xfxf (i = 2,3) as shown in Fig. [121 The left plot shows the 
results as a function of Mh±, whereas in the right plot we present the decay widths as a 
function of (pMi in S4. We show separately the results for the ^2 and h^ decay widths. In 
the left plot of Fig. [12] the hrst “apparently single” dip in the h 2 decay (upper lines) is in 
reality coming from the thresholds TN6 and TN7, see Tab. [5l The second (large) dip is the 
threshold TNIO. The last “apparently single” dip is in reality coming from the thresholds 
TNll and TN12. The “step” (anomalous threshold) at Mh± ~ 1310 GeV could be traced 
back to the G-functions Go,i, 2 (m.^±,m^±,m^,m? with s ^ s'. Away from the 

production threshold relative corrections of ~ —3% are found in S4 (see Tab. [3]) for the /i 2 
decay. The loop corrections increase with increasing Mii± and reach ~ —10% in S5. In 
case of the hs decay the dips are the same as for /i 2 and the relative corrections are only 
~ —2% in S4 (see Tab. [3]). The two mass crossings MGl and MG2 are clearly visible at 
Mii± 1006 GeV and Mh± 1532 GeV as described above, where ^2 and hs change their 
role. Between MGl and MG2 we hnd F(h 2 —)■ xfx^) > —)■ xfxf): outside it is vice 

versa, as can be clearly observed in the left plot of Fig. [T21 The suppression of the CP-even 
decay (lower lines) vs. the CP-odd decay (upper lines) is clearly visible, where at threshold 
the behavior follows Eqs. m, dini). After the threshold the decays grow roughly linear with 
the Higgs boson masses. 

We now turn to the phase dependence of the decay width shown in S4, i.e. for Mh± = 

our case the Z-factor matrix is given by Zij = ZHiggs[j, j], see Ref. [IS] (and Ref. [H]), which 
contributes at tree-level. Furthermore Z is calculated by FeynHiggs which uses mb(mb) and tree-level 
sfermion masses instead of the shifted masses, causing a slight displacement in the threshold position. 
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r/GeV hi^xfx^ r/GeV hi ^ x^xf 



Figure 12: T{hi —)■ Tree-level and full one-loop corrected partial decay widths are 

shown. The left plot shows the partial decay width with Mh± varied. The right plot shows the 
complex phase varied with parameters chosen according to S4 (see Tab. [3]). 


1200 GeV, where the right plot in Fig. [T2] shows the dependence of T{hi —)■ xfx^) V^Mi- 
Since Mi does not appear in the chargino mass matrix, the effects of varying its phase enter 
only via loop corrections and are extremely small. The relative corrections in S4 are the 
same as in the left plot. 

The results for T{hi —)■ ), as shown in the left plot of Fig. [TTlare smaller by roughly 

a factor of 2 w.r.t. F(hj —)■ xfxf)^ largely related to the kinematic suppression. At Mh± = 
1400 GeV the full one-loop corrections to the ^2 decay reach ~ -|-4%. For the decay of the 
hs aX Mh± = 1400 GeV we hnd full corrections at the level of less than -|-1%. As in the 
upper left plot one can observe the MG2 with an “interchange” of and h^. The same 
suppression of the CP-even vs. the CP-odd decay, as in Fig. [THis clearly visible. 

In the right plot of Fig. [T3]we show the complex phase varied at Mh± = 1400 GeV. 
The variation with ipMi is extremely small (for the same reasons as explained above), there¬ 
fore the full relative corrections in S5 are the same as in the left plot, see above. 

The results for the “mixed” decay, F(hj —?> Xi^x^), are shown in Fig. [141 where in the 
left (right) plot we show the dependence on Mh± (<^Mi)- In the left plot the hrst dip 
in the ^2 decay (lower lines) is the threshold TNIO, see Tab. O The remaining dip (at 
Mjh± 1268 GeV) is caused by the two thresholds TNll and TN12. At Mjh± = 1200 GeV 
the full one-loop corrections to the h 2 decay reach ~ -|-20%. Now we turn to the corresponing 
hs decay. The hrst dip (hardly visible in the upper lines) is the threshold TN8. The second 
dip is the same as for the ^2 decay, see above. For the decay of at Mh± = 1200 GeV we 
hnd full corrections at the level of -|-10%. As in Fig. [T2]one can observe the MG2 with an 
“interchange” of ^2 and h^- 
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Table 5: Thresholds in neutral Higgs boson decays. 


TNI 

Mh± 

TN2 

Mh± 

TN3 

Mh± 

TN4 

Mh± 

TN5 

Mh± 

TN 6 

Mh± 

TN7 

Mh± 

TN 8 

Mh± 

TN9 

Mh± 

TNIO: Mh± 

TNll: Mh± 

TN12: Mh± 

TN13: Mh± 


805 GeV 
^ 948 GeV 
954 GeV 
1092 GeV 
^ 1107 GeV 
^ 1108 GeV 
^ 1112 GeV 
1138 GeV 
1168 GeV 
^ 1171 GeV 
1268 GeV 
^ 1268 GeV 
1545 GeV 


m^o + m^o 
m-±+m^± 
m^o + m^o 

+ m-b^ 
m^± + m^± 

Xl X2 

m-± + m^± 

Xl X2 

m^o + m^o 

X2 X4 

Tft-zO Tft-^O 
X3 A4 

mi^ + mi^ 
mil + 
m-± + m-± 

X2 X2 

mi, + mi. 


mh 2 ~ 799 GeV 
m /13 ~ 945 GeV 
m /13 ~ 951 GeV 
m/i 2 1086 GeV 
mhs ~ 1105 GeV 
m/i 2 ~ 1105 GeV 
m/i 2 ~ 1108 GeV 
mh^ ~ 1135 GeV 
m/i 3 = 1165 GeV 
m/j 2 = 1165 GeV 
m/i 2 ~ 1264 GeV 
m /,2 ~ 1265 GeV 
m /,2 = 1542 GeV 


In the right plot of Fig. [14] one can see that the variation with ifMi is again very small 
with tiny CP-asymmetries and the same corrections as in the left plot (for the same reasons 
as explained above). 

Overall, for the neutral Higgs decays to a chargino pair we observe, again as expected, an 
increasing decay width oc Mh±, as m /,23 increase nearly linearly with our input parameter 
Mh±. The full one-loop corrections reach a level of 10% for decay widths being of 0{1 GeV), 
and they can reach up to 20% in the “mixed” decay mode. The variation with (pMi is found 
to be negligible, as expected, since Mi enters only via the loop corrections. 
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r/GeV r/GeV 



Figure 13: T{hi —)■ xjx^)- Tree-level and full one-loop corrected partial decay widths are 
shown. The left plot shows the partial decay width with Mh± varied. The right plot shows the 
complex phase ipMi varied with parameters chosen according to S5 (see Tab. [3]). 


r/GeV hi^xfxf r/GeV /j. 



Figure 14: T{hi —)■ Tree-level and full one-loop corrected partial decay widths are 

shown. The left plot shows the partial decay width with Mh± varied. The right plot shows the 
complex phase varied with parameters chosen according to S4 (see Tab. [3]). 
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hi decays into neutralinos 

Next we consider hi decays into neutralinos with equal indices. First, we present the decay 
hi —)■ XiXi ill Fig- [El Bounds on m^o often assume an underlying GUT based on a simple 
Lie group, leading to Mi = 5/3(sw/cw)^M2. If the latter assumption is dropped, hardly any 
direct bound on m^o can be placed [SB]. Therefore, we also treat Mi as an independent 
parameter. The decay hi —)■ XiXi constitutes an invisible decay of the light Higgs boson, 
where bounds based on a combination of LHC and Tevatron data yield an upper bound 
on an invisible branching ratio of ~ 40% [HI]. Since the Xi constitutes a perfect Dark 
Matter candidate in the MSSM [35], in order to determine the Dark Matter properties a 
precision measurement of this process at the LHC or a future e’''e“ collider will be necessary. 
Consequently, a precise prediction of r(hi —)■ x^xT) particular interest. 

In the upper left plot of Fig. [15] we show the results for Mh± varied in SI, but with 
I Ml I = 50 GeV as the base scenario. The full loop corrections are ~ +25% at Mh± = 
700 GeV in the upper left plot. ipMi is varied in the upper right plot. One can observe a 
strong dependence of the decay width on which can change by a factor of 8. The largest 
loop corrections are found as ~ +31% for = 72°, 288° and ~ +59% at ipM^ = 180°. In 
the lower left plot of Fig. [15] we show the decay width with Mi varied. Close to Mi = 0 the 
lightest neutralino becomes massless. For not too small values a decay width of 10-^ GeV 
can be observed, going to zero at the kinematic threshold. The one-loop corrections reach 
up to ~ +30% at Ml = 20 GeV. Finally, in the lower right plot |+| is varied, and the decay 
width drops down to ~ 10“^ GeV for /i > 600 GeV and with corrections between ~ +6% 
and ~ +28%. 

We now turn to the decays of the heavy neutral Higgs bosons. In Fig. [16] we present the 
results for the decays hi —>■ XiXi (* = 2,3). The dependence on Mh± is shown in the left 
plot, whereas the dependence on for Mh± = 700 GeV is given in the right plot. We 
start with F(hj XiXi) i^i fhe left plot. The hrst dip (lower lines) in the h 2 decay is the 
threshold TNI, see Tab. [5] The second dip (hardly visible in the upper lines) is the threshold 
TN40 The third “apparently single” dip is (again) in reality coming from the thresholds 
TN6 and TN7. The fourth (large) dip is the threshold TNIO. The last “apparently single” 
dip is in reality coming from the thresholds TNll and TN12. The full loop corrections are 
~ +11% at Mh± = 700 GeV. Also shown in this plot is the decay h^ XiXi- The hrst dip 
(upper lines) is in reality coming from the thresholds TN2 and TN3. The second dip (lower 
lines) in the h^ decay is the threshold TN5. The third dip (lower lines) is the threshold 
TN8 and the last dip is the threshold TN9. The full relative corrections reach ~ +10% at 
Mh± = 700 GeV. The suppression of the CP-even decay (lower lines, going with Eq. flTOl) at 
threshold, and then roughly linear with m/ 12 ) f^e CP-odd decay (upper lines, going with 
Eq. ([9]) at threshold, and then roughly linear with m/^g) is again clearly visible. 

In the right plot of Fig. [16] we show the h 2 decay with the complex phase cpMi varied at 
Mfj± = 700 GeV. The variation with ipMi is found to be very large, changing the decay width 

It should be noted that this threshold enter into the tree-level only via the Z matrix contribution. 
These effects propagate also into the loop corrections via 2 RefAI *ree-Mi.ioop}- Furthermore Z is calculated 
by FeynHiggs which uses mbijnb) and tree-level sfermion masses instead of the shifted masses, causing a 
slight displacement in the threshold position. 
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r X lO-VGeV hi ^ x?X? r x 10-VGeV hi x°iX°i 



r X lO-^GeV hi ^ x?X? r X lO-^GeV hi ^ x?X? 



Figure 15: r(/ii XiXi)- Tree-level and full one-loop corrected partial decay widths are 
shown with parameters chosen according to Si (see Tab. [3]) but here with \Mi\ = 50 GeV. The 
upper left plot shows the partial decay width with Mh± varied. The upper right plot shows the 
complex phase ^pmi varied. The lower left (right) plot shows Mi (|yu|) varied. 


by up to a factor of 5 where the full relative corrections are up to ~ -|-20% at ipMi = 180° 
for SI. The decay with the complex phase lpmi shows also a very large variation at 
Mh± = 700 GeV and the loop corrections reach up to ~ -t-19% at tpMi = 180°. 

In Fig. [T7]we show the decays hi —)■ analogy to Fig. [161 The same behavior 

of h 2 and concerning MGl and MG2, as well as the CP-even and CP-odd decay can be 
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observed. The dips (some are hardly visible) are the same as already described in Fig. [16] 
beginning at Mh± 1092 GeV, see above. The full relative corrections for the ^2 {h^) decay 
are ~ —18% (~ +10%) at Mi^± = 1200 GeV, i.e. S4. 

In the right plot of Fig. [T7|we show the variation of F(h 2 —)■ X 2 X 2 ) with ipMi at Mh± = 
1200 GeV. Here the loop corrections can vary between ~ —18% for 99^1 = 0°, 360° and 
~ —13% at ifMi = 180°. The decay with 99^1 varied reach ~ +11% for 99^1 ~ 90°, 270° 
in S4. 

Next, in Fig. [TS] we present the decays hi —?> + 3 + 3 , in full analogy to Fig. [ini The 
same behavior of h 2 and concerning MGl and MG2 and the CP-even/-odd decay can 
be observed. The dips (some are hardly visible) are again the same as described in Fig. [16] 
beginning at Mh± ~ 1092 GeV, see above. The “knee” at Mh± ~ 1545 GeV (red line) 
is the threshold TN13 (see Tab. [5]) in the G-functions Go,i,2(’^|o, 

The full one-loop corrections for the /i 2 decay are ~ +172% at Mh± = 1200 GeV. This 
strange behavior is a numerical effect caused by an interplay of anomalous thresholds in 
Co,i, 2 (’^|o,^^^ 2 ’('® 7^ with the effects induced by the Higgs mixing 

matrix Z. This effect is absent in the decay of the ^ 3 , where we hnd the full relative correc¬ 
tions at the level of ~ +25% for Mh± = 1200 GeV. 

In the right plot of Fig. [18] we show the variation of F(h 2 —?■ + 3 + 3 ) with at Mh± = 
1200 GeV (i.e. at an “extreme” point for the ^2 decay). Here (for the same reasons as in 
the left plot) the loop corrections reach ~ + 111 % at tpui = 180°. Also in the right plot of 
Fig. [18] we show F(h 3 —)■ X 3 X 3 ) with ipMi varied in S4. Here the loop corrections can reach 
~ +26% at if Mr = 90°, 270°. 

In Fig. [19] we present the decays hi —)■ + 4 + 4 , again in full analogy to Fig. [16] The same 
behavior of h 2 and h^ concerning MG2 and the CP-even/-odd decay can be observed. The 
full relative corrections for the ^2 decay are ~ +4% at Mu± = 1400 GeV, i.e. S5, while the 
/i 3 decay shows relative corrections less than + 1 % at Mh± = 1400 GeV. 

In the right plot of Fig. [19] we show F(h 2,3 ^ X 4 X 4 ) at Mh± = 1400 GeV. For both 
decays the variation of fMi is very small, i.e. the loop corrections reach the same values (in 
S5) as in the left plot. 

We now turn to the neutral Higgs decays to neutralinos with different indices. In this 
case, contrary to the decay into identical charginos, the CP-asymmetries are also zero, due 
to the Majorana nature of the neutralinos. 

In Fig. [20] we present the decay hi —)> XiX 2 - la the left plot we show the results as a 
function of Mjj±. The dips are (again) the same as already described in Fig. [16] beginning 
at Mff± 948 GeV, see above. The full relative corrections for the /i 2 decay are ~ +9% at 
Mff± = 900 GeV (i.e. S2). The full one-loop corrections for the h .3 decay at Mh± = 900 GeV 
reach +7%. 

In the right plot of Fig. [20] we show the hi decay with the complex phase fMi varied at 
Mh± = 900 GeV. The variation with fMr is found to be very large and the loop corrections 
vary between ~ +9% for fMi ~ 0° and ~ + 6 % at fur = 90°. We also show the hy, 
decay in the right plot of Fig. [20] with fUr varied at Mh± = 900 GeV. The variation 
with ifMi is yet larger than in the /i 2 case. The full corrections can reach ~ +10% at 
ifMr = 180°. The general behavior can be understood as follows. For fMr = 0° one hnds 
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CV{h 2 x\')<^ = ~CV{h^x\')^) = +1) leading to the above discussed suppression of the h 2 
decay, see Eqs. ([H]), ffTU]) . Going to = 180° changes the CP-nature of the Xi: leading to 
C'P{h 2 XiX 2 ) — = “1 corresponding suppression of the decay. 

In Fig. |2T]we present the decay hi —)> x?X 3 . It should be noted that the decay hi —)> y'j’xl] 
looks quite similar to Fig. [20] but with an interchange of h 2 with h^. In the left plot we show 
the results as a function of Mh±- The dips are (again) the same as already described in 
Fig. dSl beginning at Mh± ^ 948 GeV, see above. The full one-loop corrections for the h 2 
decay reach ~ +7% at Mh± = 900 GeV (i.e. S2). The relative corrections for the h^ decay 
at Mu± = 900 GeV are -1-14%. In comparison with Fig. [201 one can observe an “inversion” 
of the relative size of the decays widths of the h 2 and the h^ (green/red lines vs. blue/purple 
lines). This “inversion” is due to the fact that CV{x\xV) = ~C-V{x\xV}- 

In the right plot of Fig. [21] we show the hi decay with the complex phase ipMi varied at 
Mh± = 900 GeV. This variation is (again) found to be very large, as can be seen in the 
right plot. The loop corrections for the h .2 decay can reach ~ -|-9% at = 180°. In the 
right plot of Fig. [21] we show also the h^ decay with 99^1 varied at Mh± = 900 GeV. There 
the loop corrections vary between ~ -|-14% for ~ 0° and ~ -1-7% at (fMi = 90°. Again 
the “inversion” (as in the left plot) can be observed. 

In Fig. [22] we present the results for T{hi —)■ XixT) a function of Mh± in the left plot. 
The tree-level decay width is accidently very small for the parameter set chosen, see Tab. [3] 
Because of this smallness, the relative size of the one-loop correction becomes larger then the 
tree-level, and can even turn negative. Therefore, in this case we added |A4i_ioopP to the full 
one-loop result to obtain a positive decay width. The dips are (again) the same as already 
described in Fig. [16] beginning at Mh± ^ 948 GeV, see above. The anomalous thresholds 
(“steps” in the red line) could be traced back to the G-functions at 


Mh± ~ 1020 GeV 
Mh± ~ 1026 GeV 
Mh± ~ 1031 GeV 
Mh± ~ 1035 GeV 
Mh± ~ 1182 GeV 


Co, 1,2 , m|o, m|o, m|o, m|o, M|), 

Co, 1,2 , m|o, m|o, m^±, ±, M ^), 

{ 2 2 2 22 2\ 
<-0,1,2 , m-0, : 

Co, 1,2 , m|o, m|o, m|o, m|o, J , 

Co, 1,2 (m- 0 , , m- 0 , ), 


(in this order). The full relative corrections for the h .2 (h^) decay are ~ —77% (~ —93%) at 
Mh± = 1200 GeV (i.e. S4). 

In the right plot of Fig. [22] we show T{hi XiXT) with the complex phase (pMi varied at 
Mh± = 1200 GeV. Here (for the same reasons as in the left plot) the loop corrections can 
be larger than the tree-level (and for consistency with the left plot we also add |A4i_ioopP 
here). The loop corrections for the ^2 decay vary between ~ —77% at ifMi ~ 0°,360° and 
~ —44% at ifMi = 180°. The loop corrections for the h^ decay vary between ~ —93% at 
(fMi ~ 0°, 360° and ~ —78% at (fMi = 180°. 

In Fig. [23] we present the decay hi X^Xs- fhe left plot we show the results as a 
function of Mh± . The dips (some are hardly visible) are again the same as already described 
in Fig. [16] beginning at Mui ^ 1092 GeV, see above. The full relative corrections for the ^2 
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decay reach up to ~ +59% at Mh± = 1200 GeV (i.e. S4). The loop corrections for the 
decay at Mh± = 1200 GeV are ~ —14%. 

In the right plot of Fig. [22] we show the hi decay with the complex phase varied 
at Mh± = 1200 GeV. For the /i 2 decay the variation of ipM-^ is very small, i.e. the loop 
corrections reach (in S4) the same values (~ +59%) for all 9 ?Mi- The variation of ipMi in the 

decay is also small with corrections at the level of ~ —14%. 

In Fig. [21] we show the decay hi —)■ x^xl- In the left plot we show the results as a function 
of Mh±. The dips are (again) the same as already described in Fig. [12] beginning at Mh± ~ 
1171 GeV, see above. The full relative corrections for the /i 2 decay reach up to ~ +15% at 
Mh± = 1200 GeV, i.e. S4. The loop corrections for the hs decay at Mh± = 1200 GeV are 
~ +9% 

In the right plot of Fig. [2l]the hi decay is shown with the complex phase (fMi varied in 
S4. For both decays the variation of (fMi is very small, as expected, since X 2 and x1 are 
determined largely by M 2 and /i in the neutralino mass matrix (for the parameters chosen as 
in Tab.[3|). The loop corrections for the ^2 (hs) decay reach ~ +16% (~ +9%) at <pmi = 180°. 

The final decays involving neutralinos are shown in Fig. [251 The results as a function of 
Mh± are given in the left plot. The dips are (again) the same as already described in Fig. [16] 
beginning at Mh± ^ 1171 GeV, see above. The full relative corrections are only ~ +3% at 
Mh± = 1200 GeV (i.e. S4). The full relative corrections at Mh± = 1200 GeV reach +6%. 
In comparison with Fig. [24] one can observe (again) an “inversion” of the relative size of the 
decays widths of the ^2 and the h^, due to the fact that CV{x 2 '^ 4 ) — ~^'P{X 3 XT>- 

In the right plot of Fig. [25] we show the hi decay with the complex phase (fMi varied in S4. 
For both decays the variation of ifMi is again very small, since both neutralinos are largely 
determined by jj, (for the parameters chosen as in Tab. [3|). The full one-loop corrections are 
the same as for the left plot. 

Overall, for the nentral Higgs decays to a nentralino pair we observed, again as expected, 
an increasing decay width oc The full one-loop corrections reach a level of 10 — 20% 

for decay widths of 0(1 GeV). The variation with (pMi is fonnd largest in cases where the 
CP-nature of the decay depends strongly on the phase, there then changes by a factor of 5 
or more can be observed. 


Exceptions are the hi XiXi decay (see the upper left plot of Fig. [15]), since rrihi depends only very 
weakly on M^±. The next exception are the corrections to the ^2 decay in the left plot of Fig. [TS] (red line), 
due to an accidental interplay of anomalous thresholds with the effects induced by the Higgs mixing matrix 
Z. The final exception can be observed in the corrections to the ^2 decay in the left plot of Fig. [22] (red 
line), because of the additional 2-loop corrections |Mi-ioopP (see the discussion of Fig. [22] above). 
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r/GeV h, ^ x?X? r/GeV h, ^ 



Figure 16: r(hi —>■ XiX?)- Tree-level and full one-loop corrected partial decay widths are shown. 
The left plot shows the partial decay width with Mff± varied. The right plot shows the complex 
phase (fiMi varied with parameters chosen according to Si (see Tab. [3]). 


r/GeV hi X2X2 r/GeV hi X2X2 



Mh± 




Figure 17: V{hi —)■ X 2 X 2 )- Tree-level and full one-loop corrected partial decay widths are shown. 
The left plot shows the partial decay width with Mh± varied. The right plot shows the complex 
phase ^pmi varied with parameters chosen according to S4 (see Tab. [3]). 
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Figure 18: T{hi XsXs)- Tree-level and full one-loop corrected partial decay widths are shown. 
The left plot shows the partial decay width with Mh± varied. The right plot shows the complex 
phase ifMi varied with parameters chosen according to S4 (see Tab. [3]). 


r/GeV ft, ^0^0 r/GeV ft, 





Figure 19: V{hi —)■ X 4 X 4 )- Tree-level and full one-loop corrected partial decay widths are shown. 
The left plot shows the partial decay width with Mh± varied. The right plot shows the complex 
phase ^pmi varied with parameters chosen according to S5 (see Tab. [3]). 
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r/GeV hi ^ x°iX°2 r/GeV ft, ^ 



Mh± V^Mi 


Figure 20: T{hi XiX^)- Tree-level and full one-loop corrected partial decay widths are shown. 
The left plot shows the partial decay width with Mh± varied. The right plot shows the complex 
phase ifMi varied with parameters chosen according to S2 (see Tab. [3]). 


r/GeV hi x°x° r/GeV hi x?X3 



Mh± 




Figure 21: V{hi —)■ XiXs)- Tree-level and full one-loop corrected partial decay widths are shown. 
The left plot shows the partial decay width with Mh± varied. The right plot shows the complex 
phase ^pmi varied with parameters chosen according to S2 (see Tab. [3]). 
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Mh± V^Mi 


Figure 22: T{hi XiXa)- Tree-level and full one-loop corrected partial decay widths are shown. 
The left plot shows the partial decay width with Mh± varied. The right plot shows the complex 
phase ifMi varied with parameters chosen according to S4 (see Tab. [3]). 


r/GeV hi -5- X2X3 r/GeV hi X2X3 



Mjj± ^Mi 


Figure 23: T{hi —)■ x^xD- Tree-level and full one-loop corrected partial decay widths are shown. 
The left plot shows the partial decay width with Mh± varied. The right plot shows the complex 
phase ^pmi varied with parameters chosen according to S4 (see Tab. [3]). 
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r/GeV hi ^ x°xO r/GeV h^ ^ x°2X°4 



Mh± V^Mi 


Figure 24: r(/ij ^ X2X4)- Tree-level and full one-loop corrected partial decay widths are shown. 
The left plot shows the partial decay width with Mh± varied. The right plot shows the complex 
phase varied with parameters chosen according to S4 (see Tab. [3]). 


r/GeV h^xlxl r/GeV hi^^xl 



Figure 25: V{hi —)■ XsX^)- Tree-level and full one-loop corrected partial decay widths are shown. 
The left plot shows the partial decay width with Mh± varied. The right plot shows the complex 
phase ^pmi varied with parameters chosen according to S4 (see Tab. [3]). 


33 





































5 Conclusions 


We evaluated all partial decay widths corresponding to a two-body decay of the MSSM 
Higgs bosons to charginos and neutralinos, allowing for complex parameters. In the case of 
a discovery of additional Higgs bosons a subsequent precision measurement of their properties 
will be crucial determine their nature and the underlying (SUSY) parameters. In order to 
yield a sufficient accuracy, one-loop corrections to the various Higgs-boson decay modes have 
to be considered. In this work we take another step in the direction of completion of the 
calculation of all two-body decays at the one-loop level in the cMSSM in this stable and 
reliable renormalization scheme: we calculated all two-body decay modes of the Higgs bosons 
to charginos and neutralinos in the cMSSM. 

The decay modes are given in Eqs. ([I]) - (I2])- The evaluation is based on a full one-loop 
calculation of all decay channels, also including hard and soft QED radiation. We restricted 
ourselves to a version of our renormalization scheme which is valid for \Mi\ < IM 2 I, |/i| and 
M 2 7 ^ fi (where Mi and M 2 denote the soft SUSY-breaking parameter of the U{1) and SU{2) 
gauginos, and p is the Higgs mixing parameter) to simplify the analysis, even though our set¬ 
up allows to switch to other parameter regions, possibly implying a different renormalization, 
see the discussion in Refs. pHHTH] . 

We hrst reviewed the relevant sectors including some details on the one-loop renormal¬ 
ization procedure of the cMSSM, which are relevant for our calculation. In most cases we 
follow Ref. [l5]. However, in the scalar fermion sector, where we differ from Ref. [15] the 
relevant details are indicated. We have discussed the calculation of the one-loop diagrams, 
the treatment of UV and IR divergences that are canceled by the inclusion of (hard and 
soft) QED radiation. We have checked our result against the literature, and in most cases 
we have found good agreement, once our set-up was changed to the one used in the existing 
analyses. 

While the analytical calculation has been performed for all decay modes to charginos 
and neutralinos, in the numerical analysis we mostly concentrated on the decays of the 
heavy Higgs bosons, with hi —)■ XiXi being the only channel analyzed for the light neutral 
Higgs boson. For the analysis we have chosen a parameter set that allows simultaneously a 
maximum number of two-body chargino/neutralino decay modes. In the analysis either the 
charged Higgs boson mass or the phase of Mi has been varied. For Mh± we investigated 
an interval starting at Mh± = 600 GeV up to Mu± = 1.6 TeV, which roughly coincides 
with the reach of the LHC for high-luminosity running as well as an e’''e“ collider with a 
center-of-mass energy up to ^/s ~ 3 TeV. 

In our numerical scenarios we compared the tree-level partial decay widths with the full 
one-loop corrected partial decay widths. We concentrated on the analysis of the decay widths 
themselves, since the size of the corresponding branching ratios (and thus the size of their 
one-loop effects) is highly parameter dependent. 

We found sizable corrections of ~ 10% in many channels, sometimes going up to ~ 20%. 
Even larger corrections are only found in cases where the tree-level result is accidentally 
small and thus the decay likely not observable. Corrections at the 10 — 20% level have 
also been found for the decay hi —)■ XiXi, which could constitute an important channel 
for the determination of the Dark Matter properties in the cMSSM. Consequently, the full 
one-loop corrections should be taken into account for the interpretation of the searches for 
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charginos/neutralinos as well as for any future precision analyses of those decays. 

The tree-level results, but also the size of the full one-loop corrections often depend 
strongly on the complex phase analyzed, The one-loop contributions can vary by a 

factor of ~ 2 as a function of the complex phase. Neglecting the phase dependence could 
lead to a wrong impression of the relative size of the various decay widths. 

In the cases where a decay and its complex conjugate hnal state are possible, i.e. the 
charged Higgs decays we have evaluated both decay widths independently. The asymmetries, 
as a byproduct of our calculation, turn out to be relatively small, at the level of a few per¬ 
cent. 

The numerical results we have shown are, of course, dependent on the choice of the SUSY 
parameters. Nevertheless, they give an idea of the relevance of the full one-loop corrections. 
Decay channels (and their respective one-loop corrections) that may look unobservable due 
to the smallness of their decay width in our numerical examples could become important if 
other channels are kinematically forbidden. Following our analysis it is evident that the full 
one-loop corrections are mandatory for a precise prediction of the various branching ratios. 
The full one-loop corrections should be taken into account in any precise determination 
of (SUSY) parameters from the decay of (heavy) MSSM Higgs bosons. It is planned to 
implement the evaluation of the branching ratios of the (heavy) Higgs bosons into the Fortran 
code FeynHiggs, together with an automated choice of the renormalization scheme valid for 
the full cMSSM parameter space. 
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